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The semiclassical model of relativistic Coulomb excitation is studied in situations 
in which the impact parameter is small enough so that projectile and target charge 
distributions overlap. The electromagnetic effects of this overlap are shown to be 
small. Realistic nucleon-nucleon reaction cross-sections, and realistic nuclear radial 
charge and matter distributions are used to determine a formula for the lower impact 
parameter limit to be used in the calculation of the Coulomb excitation cross-section. 
A wide selection of projectile-target pairs is explored, in the bombarding energy range 
of 1 GeV to 5 GeV per nucleon. 



The semi-classical approach to relativistic Coulomb excitation (RCE) assumes that the 
Coulomb impulse suffered by the projectile nucleus as it passes the target nucleus is small 
compared to the linear momentum of the projectile. Then the deflection of the projectile will 
be small, and its trajectory can be approximated by a straight line, with impact parameter 
b. If b is large enough so that the projectile and target do not experience each others 
nuclear forces, the only processes that can occur are electromagnetic. For example, as the 
projectile moves along its trajectory, the electromagnetic fields it produces can induce induce 
transitions in the target [l^. The RCE cross-section for the population of a final target state 
ip a , starting with the target ground state ip g . s ., is given by 
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where P g . s .-> a (b) is the probability of the ipg.s. —> i>a target transition when the target 
experiences the electromagnetic impuse due to the projectile following an orbit with impact 
parameter b. 

The choice of b m [ n can have a significant effect on (T g , s ,_» a calculated with Equation (jl.lj) . 
since P g . s .-> a (b) has its greatest value at b = 6 min . This is because the projectile electromag- 
netic fields at the target are greatest when the two are in closest proximity. Also, for small 
b the electromagnetic impulse at the target is more sudden, which makes it more effective 
at exciting high-energy target states, such as giant multipole excitations. One of the goals 
of this paper is to provide reliable values of 6 m i n for a wide variety of projectile-target pairs 
and a wide range of relativistic bombarding energies. 

The assumption behind Equation (jl.lj) is that there are no nuclear interactions when 
b > fc m in; but when b < 6 min the nuclear interactions are so strong that they dominate over 
electromagnetic processes, and we no longer have Coulomb excitation. We will refer to 
this latter situation using the term absorption. A more realistic discussion would involve a 
continuous transition as b increases, from complete absorption at small b, to zero absorption 
at large b. If X(b) is the probability that a projectile traverses its orbit with no nuclear 
interaction with the target, the absorption probability would be 1 — X(b), and the RCE 
cross-section would be 



The function X(b) can be estimated by folding the nuclear density distributions of the 
projectile and target with the nucleon-nucleon interaction cross-section. 

For the small-6 range of the integral in Equation (jl.2jl there are situations in which the 
tails of the projectile and target overlap. When this happens, the electromagnetic interaction 
between the projectile and target is more complicated than the situation of Equation (jl.lj) . 
where all the projectile charge density is assumed to be outside the target. In this connection, 
we explore the 6-range of the integral in Equation (jl.2j) where X(b) is making the transition 
from a small value (almost complete absorption) to near unity (almost no absorption). 

More precisely, in Sec. II we consider the nucleon-nucleon interactions, and an explicit 
form will be deduced for the function X(b). In Sec. Ill we will investigate the effect on 
Pg.s.^a{b) of the overlapping of the projectile and target charge densities. We will then 
discuss in Sec. IV several prescriptions for the parameter 6 min of Equation (jl.lj) which will 
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lead to the same calculated cross-section as the more accurate Equation (jl.2j) . 

II. NUCLEAR INTERACTION BETWEEN PROJECTILE AND TARGET 

We discuss the nuclear interactions between the projectile and target nuclei in terms of 
the optical limit of the Glauber model yj, in which each projectile nucleon is treated as a 
grey disk, characterized by a profile function 7(3). In Figure 1, the projectile and target 
centers are labelled O p and T , respectively. The edge view of a disk is shown, representing 
a projectile nucleon located at rp relative to Op. Points in the plane of the disk are located 
by s relative to the center of the disk (the center of the projectile nucleon). If the projectile 
moves a distance dz, then the disk sweeps out an effective reaction volume of 

dz J d 2 s j(s), (2.1) 

where ■y(s) is the profile function, related to the nucleon-nucleon scattering amplitude /nn(q) 
by a two-dimensional Fourier transform 

7(8) = J d 2 q e-^ s / NN (q). (2.2) 

If n T (r T ) is the number density of target nucleons, then the probability that a target nucleon 
is in the effective reaction volume given by Equation (j2.1j) is 



dz J d s 7(5) n T (r T ), 

where 



r T = b + zz + r p + s. (2.3) 

This refers to a single projectile nucleon. If n p (r p ) is the number density of projectile 
nucleons, the total probability of a nuclear interaction when the projectile moves a distance 
dz is 

i/ft f M.,)/Ar(.)^,)^/H. 

Therefore the probability P(b) that no nuclear interactions occur when z goes from —00 to 
00 is 

X(b) = e~f-°° dz = e~ dz f d3rp n p( rp) 1 d2s ^ n T( r T) (2.4) 



Equation ()2.3|) implies that, for fixed r p and b, 

d 2 s dz = d 2 (r T )± dz T = d 3 r T 
s = (r T )± - (r p )_L - b 

Thus Equation (|2.4|) can be re-written as 

X(b) = e~^ d3r T d3r T "p^p^t^t) 7(l(r T )_L-(r p )x-b|) (2.5) 

This is the X(b) needed in Equation f)1.2j) . 

To evaluate the integral in the exponential of Equation (J2.5jl . we express the number 
densities in terms of their Fourier transforms. For spherically symmetric nuclei we define 

f An f°° 

n p (q) = / d 3 r p n p (r p )e iq ' rp = — / r sm(qr)n p (r)dr = rip(q) (2.6a) 

J q Jo 

f Ait f°° 

n T (q) = / d 3 r T ra T (r T )e iq ' rT = — / r sm(qr)n T (r)dr = n T (q) (2.6b) 
J Q Jo 



With the help of Equation (|2.2jl . we can write 

J d 3 r p d 3 r T n p (r p )n T (r T ) 7(|(r T )j. - (r p ) ± - b|) 

d 2 q ± e^ b n p (q ± )n T (q ± ) f NN (q ± ) 



1 



2nik NN 



1 f°° 

— / q±dq ± J (q L b) n p (q L )n T (q L ) f NN (q±) (2.7) 

lK NN JO 



We have so far not distinguished between neutrons and protons. Since scattering ampli- 
tudes for identical and non-identical nucleons are different, the integral in Equation ()2.5|) 
should be replaced by 

J d 3 r p d 3 r T [n pz (r p )n TZ (r T ) + n P N (r p )n TN (r T )] 7zz(|(r T )± - (r p ) ± - b|) 

+ J d 3 r p d 3 r T [n p z (r p )n TN (r T ) + n PN (r p )n TZ (r T )] 7z N (|(r T )± - (r p ) ± - b|) 

Here the subscripts Z and N distinguish between protons and neutrons, whether occurring 
in the projectile or target. Similarly, Equation ()2.2|) should be replaced by 

1 f ,2. 
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If the center-of-mass energy of the projectile and target nucleons is in the GeV region, as 
in relativistic Coulomb excitation, the spin-averaged nucleon-nucleon scattering amplitudes 
can be parameterized by (see, e.g., ref.[2|) 



e -Bzz(hq±) 2 



ik Z z 47r 

/zN(g±) _ ^ZN e -B w {hq x )* 
ik Z N 47T 

The numerical values of azz, -S zz ,o"zn and £>zn, which are functions of energy, were taken 
from the article by Igo Q|. They are given in Table 1. 

The nucleon number densities required in Equations ()2.6al2.6b|) were taken from the 
IAEA compilation [3 : ]. They are the results of Hartree-Fock-Bogoliubov calculations, whose 
parameters are fitted to measured nuclear masses. In a comparison with 523 nuclei, these 
number densities agreed with measured radii within an rms error of .028 fm. 

Bertulani et al |4j have used a similar approach to represent the effect of nuclear absorption 
in grazing collisions. Their analysis assumes that the profile function 7(5) has zero extension, 
equivalent to the assumption that the nucleon-nucleon scattering amplitude is q independent. 
They do not distinguish between the interaction of identical and non-identical nucleons. In 
different contexts, a similar approach has been developed by Benesh, Cook and Vary P and 
by Kox et al 



III. THE ELECTROMAGNETIC TRANSITION MATRIX ELEMENTS WITH 
OVERLAPPING CHARGE DISTRIBUTIONS 

In the previous Sections, we saw that at the bottom of the range of impact parameters 
there can be overlapping of the projectile and target charge distributions. We now consider 
whether this requires us to modify the standard analysis of relativistic Coulomb excitation, 
which is based on the assumption that the projectile and target charge distributions do not 
overlap. 



A. The scalar and vector potentials 

We orient our axes so that the projectile center moves with velocity vz relative to the 
target center in their common y — z plane. If (x, y, z, t) locate an event relative to the 
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projectile, and (x',y', z',t') locate the same event relative to the target, then 

o 

V V i 

x = x, y = y -b, z = i(z -vt), i = 7(t'-— z), 7 = (1 - — )s 

c c 

Suppose that the projectile charge density is static and spherically symmetric in its own 
rest frame, so that it can be written 



P P {x,y,z) = P P {r), r 



\Ac 2 + y 2 + z 2 . 



Then the scalar potential due to this charge distribution, as measured in the projectile rest 
frame, is given by 

47r r f°° 

P ( r ) = — / s 2 p p (s)ds + Aix / sp p (s)ds (3.1) 

r Js=0 J s=r 

It can be readily verified that this expression for P (r) satisfies Poisson's equation for the 
specified projectile charge density: 

V r 2 P (r) = ^ r2 ^p( r ) = -47rp P (r). 

The scalar potential measured by an observer at the target center can be determined 
from Equation ()3. 1)1 . since the scalar potential transforms as the 0-component of a Lorentz 
4-vector. This leads to 

(j)' p (x',y',z',t') = j<p p (x,y,z,t) 

47T7 



r 



pr roo 

/ s 2 p p (s)ds + 47T7 / sp p (s)ds (3.2) 

J s=0 J s=r 



with r = a/x' 2 + {%/' — b) 2 + ^y 2 (z' — vt') 2 . Note that there is no contribution in Equa- 
tion ()3.2|) from a projectile current density, as seen by a projectile-based observer. This 
would be the situation, for example, if the projectile were a doubly-closed-shell nucleus, in 
which every occupied shell-model state was matched by an occupied time-reversed state. 

It is helpful to modify Equation ()3.2|) by adding and subtracting the quantity 
— f°° s 2 p p (s)ds. This converts Equation ()3.2|) into 



4 7r 7 f 00 f°° s 2 
(j)' p (x',y',z',t') — / s 2 p p (s)ds + 47T7 / (s )p p (s)ds. (3.3a) 



s=0 



r 



Since the total projectile charge is Z p e = f^ Q 4:7rs 2 p p (s)ds, this equation can be re-written 



as 



<P' p (x', y', z', = + 47T7 (s - -) Pp (s)ds. (3.3b) 

a/x /2 + {y' - 6) 2 + 7 2 (^' - Vt') 2 J s =r r 
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The first term on the right-hand sides of Equations (j3.3al3.3bj) is the Lienard-Wiechert 
(LW) potential (see, e.g., |7(), which is the potential at (x',y',z',t') produced by a point 
charge Z p e at the projectile center, as measured by a target-based observer. The second 
term is a correction to the LW potential, needed when the projectile has charge that extends 
farther from the projectile center than the observation point (x', y', z', t'). In other words, 
this correction term is needed for points in space within the projectile charge distribution. It 
is seen that although the LW term involves only the total projectile charge Z p e, the overlap 
correction term depends upon the radial shape of the projectile charge distribution. 

The same Lorentz transformation that yields the scalar potential ()3.2j) yields the vector 
potential 

A' P (r',0 = -<(r',i')z, (3.4) 



B. Calculation of the transition charge density 

Let us consider the specific case of a direct transition between the zero- and one-phonon 
states of the target giant dipole excitation. We will generate the transition charge and 
current densities with the Goldhaber- Teller model jsj of the giant dipole resonance (GDR), 
in which the target protons and neutrons oscillate relative to each other. For example, 
phonon excitations of the 40 Ca ground state can be constructed in which spherical clusters 
of 20 protons and 20 neutrons oscillate relative to each other, with the oscillation degree of 
freedom the vector r pn , drawn from the center of the neutron cluster to the center of the 
proton cluster (see Figure 2). If the oscillations are small, the GDR potential can be expected 
to be approximately harmonic, and then the relative motion of the cluster centers would be 
governed by a harmonic oscillator wave function ). The ground state relative motion 

would be determined by ?/>o' (r pn ), and the ground state charge density would be 

Pg.s.(r') = e J d\ n [#(rj]* ^ °' Q (r pn ) G (|r' - /r pn |) , (3.5a) 

where G(s) is the number density of the protons at a distance s from the center of the 
protons (see Figure 2), and / is defined by 

f " z^TJ^ = M> (3 - 5b) 

so that fr pn is the vector connecting the target center of mass to the center of the target 
proton cluster, which we take to be an inert sphere. The transition charge density for the 



transition from the ground state to the GDR state with relative motion (n, £, m) is given by 
p g ,.^n,Lm)(r>) = e [ d\ n [rJ(rjY i#°(rj G (|r' - /r pn |) , (3.5c) 



It will be sufficient for our purposes to consider the transition from the ground state to 
the mode 

^ - ^ [vf 1 + A 1 ]. 

which is one quantum of oscillation in the y direction. At high bombarding energy, this is 
the strongest direct transition. Because this state has a simple interpretation in a Cartesian 
representation, we re-write Equations (|3.5al I3.5c|) in terms of one-dimensional harmonic 
oscillator eigenstates: 

Pg.s.(r') = e J dx pn dy pn dz pn [^oO^pn)] 2 l4>o(y P n)} 2 [4>o(z pn )] 2 

x G (x' - f x pn ) 2 + (y'-f y m f + (z' - f z pn ) 2 ^j (3.6a) 

Pg.s.^y (r') = e dx pn dy pn dz pn [ip (x pn )} [^i(?/pn)^o(l/pn)] [^o(z pn )] 

x G (j/ (x> - f x pn ) 2 + (y'-f y m f + (V - / ^n) 2 ) (3.6b) 
The one-dimensional harmonic oscillator states needed here are 

Mv P «) = e-i^ 

^(j/pa) = (^) * V 7 ^ 2/ pn e"^» = - y- ^— ^ (Z/pn) 



so that 



/2 <9 1 9 

- ^o(y P n) « ^o(Z/pn) = 7= 77 ft/>o(2/pn)]' 

^ cy P n y2u oy P n 

If we use this relation in Equation (|3.6b|) 

e r , , , r , , m2 5 



Pg.s.^y (r') = 7= / dx pn dy m dz pn [ipo(x pn )} - — [i>o(y P n)] bl>o(zpn)Y 

y2v J oy pn 

x G (J (x~> - f x pn f + {y> - f y pn ) 2 + (z'-f z pn )' 
and integrate by parts, we get 

Pg.s.^y(r') = + -= / dx pn dy pn dz pn [ip (x pn )] 2 [*Po(y P n)} 2 [^o(^ P n)] 2 



'2u 



x ^- G (yV - / x pn ) 2 + (y'-f y m ) 2 + (z - f z pn ) 2 ^j . (3.7) 
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But 







G 



yj{x'-f x pn ) 2 + {y'-f y P n) 2 + (V - f z pn ) 2 



pn 



dy' 



7 G ( \](x' - f x pn ) 2 + (y'-f y pn ) 2 + (z' - f z pn ) 2 



so that Equation ()3.7|) becomes 



Pg.s.^y 




(3.8) 



Thus the radial shape of the transition charge density is obtained from the radial derivative 
of the target ground state charge density. This result depends upon the interpretation of the 
GDR as a harmonic oscillation of proton and neutron spheres relative to each other, but it 
makes no assumption about the radial charge dependence of the proton sphere. 

The occurrence of the derivative of the ground state radial density in expressions for 
transition densities is familiar from discussions of direct inelastic scattering to modes which 
are interpreted as small shape oscillations. An exhaustive discussion of this topic can be 
found, e.g., in Ch.14 of ref. [9J. These theories are based on expansions of the nuclear shape 
in powers of the small parameter describing the oscillation. We note that our derivation of 
Equation ()3.8j) is not based on an expansion in powers of a small parameter. However, our 
use of harmonic oscillator wave functions in the derivation of Equation ()3.8|) depended upon 
our assumption that the GDR potential is harmonic, which will generally be true only if the 
deviation from equilibrium is small. 

The size parameter v used in Equations (|3.fia|) through ()3.8j) is given by v = fiu/h, where 
/i is the reduced mass for the neutron and proton spheres, and to is the GDR oscillation 
frequency. We will use the generic formula [l0( 



hu = 79 Ap 3 , 



(3.9) 



where is the target mass number. 



10 



C. The transition matrix element 

Some of the basic formulae of the semi-classical approach to RCE are presented in the 
Appendix. 

The time-dependent target transition matrix element is 



V g . s .^ y (t', b) = I dV 



,_ y (r')0 / p (r',t')-^j g .,^ y (r')-A;(r',t') 



(3.10) 



According to Equations (|3.4j) and (|3.10|) . only the z component of the transition current 
density enters into the matrix element. This is zero for the matrix element connecting the 
zero-phonon state to the state with one phonon of oscillation in the y direction. Thus, for 
this particular matrix element, we need only be concerned with the p — (ft term in ()3.10|) . 

The partition of the scalar potential into LW and overlap terms carries over to the tran- 
sition matrix element: 

V g . 8 .^ y (t',b) = V hW (t') + Kvcrlap^') (3.11a) 

f f°° S 2 

Overlap = / dV / (s ~ —) Pp (s)ds X p g . s .^ y (r') (3.11c) 

J J s=r T 

with r = ^x' 2 + (y' - b) 2 + ^ 2 {z' - vt') 2 
For the purposes of calculation, it is convenient to define a function F(r) by 

roc 

F(r) = 4tt7 / s(s - r)p p (s)ds. (3.12) 

J s=r 

Then Vlw(^) an d Kveriap(^) can be written more compactly as 

V LW (t') = [ rfV Pg - s ^ y(r V (0) (3.13a) 
J r 

V oveAap (t') = ~ J dV Pg - s ^ y(r ° F(r) (3.13b) 
with r = a/ x' 2 + (y' — b) 2 + -y 2 (z' — vt') 2 

Once p P (s), the radial form of the projectile proton density distribution, has been chosen, 
F(r) is calculated using Equation (j3.12|) . and it is then used in the numerical evaluation of 
the integrals ()3.13al3.13b|) . 
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D. Cross-sections 



Up to bombarding energies of about 5 GeV per nucleon, first-order time-dependent per- 
turbation theory accounts for almost all of the population of the state with one-phonon of 
oscillation in the y direction (see, e.g., jllj]). This leads to 



(b) 



00 dt' 



e^'V g . s ^ y (t',b) 



V g . S .^y(uJ, b) 



where V g . s .^ y (u;, b) is the "on-shell" Fourier component of V g , s ^ y (t', b) (see Appendix). Using 
Equation (jl.2|) . we can write 



a 



bdb 



s.^y 



[u,b)\ X(b) 



bdb 



V LW {U,b) + Overlap (w , 6) X (b) 



bdb [y LVf (u, 6) J X(b) + j bdb V OVCI i ap (uj, b) x [W^u, b) + V oveAap (u, 6) J X(b) 

= CLW + ^overlap- (3-14) 

Here <tlw is the cross-section that would have been calculated had the Lienard-Wiechert 
potential been used everywhere, even in the overlap region, and o" over iap is the correction that 
must be applied due to the inadequacy of Lienard-Wiechert potential. 



Application of the Winther-Alder 



12| general formula, Equation (|Alj) . to Vlw^j b) gives 



VLw (<*>,&) 



(3.15) 



h'yv 2 V v \l v ) Jo JL \c J dr' 

We can use this in Equation f)3.14|) . but Kveria P (^, b) must be evaluated numerically, as must 
the b integral in cr vcria P - The latter is simplified by the fact that the effective b range is 
finite. It is limited from below by the vanishing of X (b), and is limited from above by the 
vanishing of Kveria P (^, b), since large b implies small overlap. 

Table 2 shows a comparison of <tlw and o" over i ap for 16 and 208 Pb projectiles bombarding 
74 Ge and 202 Hg targets at 2 GeV per nucleon. It is seen that <tlw is very much larger when 
208 Pb is the projectile, and when 202 Hg is the target. This is because all the projectile and 
target charge contribute to the LW cross-section. However the amount of overlap charge 
is approximately the same in all four cases, so that the four overlap corrections are of the 
same order of magnitude. Another feature that surpresses the overlap correction in the 
heavier systems is their greater neutron/proton ratio. The extra neutrons contribute to the 
absorption, but not to the electromagnetic interaction. 
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Although <7 over i ap is a greater fraction of o"lw i n lighter systems than in heavier sstems, 
we can see from Table 2 that cr vcriap is always very small compared to olw- Thus we can 
safely ignore a" over i ap , and use <tlw alone to account for experimental data. 



IV. PROJECTILE AND TARGET DEPENDENCE OF b n 



Suppose that X (b) makes a sharp transition from to 1 as b crosses a particular value 
b min . Then Equation ()3.14j) yields 

2 



0"LW 



bdb[V LW (u,b) 



(4.1) 



Application of Equations ()A1|) through (jA4|) ) gives 



0"LW 



Tlbt 



Ulbr, 



ujb„ 



7^ 



jv 
uib m u 



ub n 



Ujbr, 



7V 



ncZpef /32 

h'jv 2 




12, / ■ (u t\ dp g .s.( r ') 
r dr ji[ —r 1 



V c 



dr' 



(4.2) 



It would be advantageous to have a prescription for 6 m i n such that the entire expression (jH.14|) 



for the cross-section would be given 
seek a formula for 6 miTn that satisfies 



3y 



16] 



the explicit formula Equation (j4.2|) . Equivalently, we 



bdb 



6=0 



V LW (cu,b) 2 X{b) 



bdb 



b=b„ 



V LW {u,bf 



Since all the b dependence of Vlw(<^, b) is contained in K\ y ) , we can write this equation 
more explicitly as 



r Ub ( Kl (^)) 2 x m = r mUA^ 

J b =o V \jvjj j b=bm . n v. V7^ 



LUb^ 



(4.3) 



We could use 6 min defined in this way in Equation ()4.2|) . thus incorporating the effect of 
nuclear interactions during grazing collisions of the projectile and target. 

Since the left-hand side of Equation (|4.3|) must be evaluated numerically, we cannot 
produce a closed formula for b m [ n . However we can numerically evaluate Equation (J4.3|) 
for a variety of projectiles, targets, and bombarding energies, and then look for regularities 
that could guide our choice of 6 m ; n in any particular situation. To carry out this program, 
we have selected as projectiles the nuclei 16 0, 40 Ca, 120 Sn, and 208 Pb, and as targets the 
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nuclei 16 0, 32 S, 52 Cr, 74 Ge, 90 Zr, 114 Cd, 138 Ba, 158 Gd, 180 Hf, and 202 Hg. The GDR oscillation 



frequencies were given by Equation ()3.9|) . For each projectile-target combination, we have 
varied the bombarding energy from 1 to 5 GeV per nucleon, in steps 1 GeV per nucleon, 
and used Equation (|4.3jl to calculate the appropriate value of b m \ n . 

In the limit of very large nucleon-nucleon interaction cross-sections, we would expect b m i n 
to be of the order of Rp + R^, or perhaps somewhat larger if the finite range of the nucleon- 
nucleon interaction is included. This would lead to an expression for 6 min proportional to 
+ Ap . However, a more realistic nucleon-nucleon interaction would allow the projectile 
and target densities to overlap slightly without nuclear interaction, which would lead to a 
smaller value of b m i n . If the projectile and/or the target have large radii, even a small 
amount of penetration implies a large overlap volume, and thus a large nuclear interaction 
probability. Thus the downward correction to an Ap 3 + Ap 3 term would be expected to 
decrease as the size of the colliding nuclei increases. The simplest way to incorporate these 
trends into a formula is to seek parameters A, /i such that 



This is the form used by Benesh et al |5J] in their analysis of the total reaction cross-section 
for colliding nuclei. We have chosen A and /x to produce the best fit, in a least-squares sense 
to the 200 (4 x 10 x 5) projectile-target-energy combinations for which we have numerically 
calculated 6 m i n using Equation (|4.3jl . The result is 



If these paramaters are used in Equation (|4.4J) . the calculated values of b min are reproduced 
with an r.m.s. error of 0.0248 fm per point. A graphical comparison is shown in Figure 
3. The continuous lines are plots of Equation (|4.4j) for our four projectiles, and the plotted 
points refer to our numerical calculations of 6 m j n for our ten targets at bombarding energies 
of 1 and 5 GeV per nucleon. It is seen that Equation (|4.4jl . with the parameters given above, 
provides an excellent representation of the 6 m i n values numerically calculated from Equation 

Another study j(| of nucleus-nucleus reaction cross-sections adopted a form equivalent to 




(4.4) 



A 



1.3115fm 



fi = 1.0509fm. 




— c 



(4.5) 
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This equation contains three adjustable parameters, A,/i and c. A least-square-deviation fit 
to our 200 values of b min yields 

A = 1.3338fm, \i = 0.1652fm, c = 1.0667fm, 

with an r.m.s error of 0.0382 fm per point. We see that the form (14. 5|) with three free 
parameters does not produce as good an overall fit as the form ()4.4j) . which has only two 
free parameters. We conclude that Equation ()4.4j) gives the most effective and economical 
representation of the dependence of b min on projectile and target mass numbers. 

The 6 m i n values given by Equations ()4.4|) or ()4.5|) yield very nearly the same CEX cross- 
sections when substituted into Equation ()4.2|) . The differences are less than 1% for all 
cases except for 1 GeV per nucleon 16 projectiles on a 16 target, where the difference is 
2%. These differences are probably small compared to the errors associated with the use of 
first-order time- dependent perturbation theory. 



V. DISCUSSION 



All the calculations presented so far have referred to /x = 1 transitions, i.e. transitions in 
which the transfer of the z component of angular momentum is ±h. These are the dominant 



transitions in RCE. Indeed, in the Weizsacker- Williams [13( approach to RCE, these are the 
only transitions considered. Nevertheless, transtions with /i ^ 1 are possible, and it is of 
some interest to know how a change in fx would affect b m \ n . This can be answered simply by 
replacing K\ in Equation ()4.3j) by K^. In the vicinity of b ~ 6 min , the argument of ( ^™ in ) 
is small, which implies that is proportional to Thus falls more sharply with 
increasing b as /x increases, and this has the consequence that b m \ n calculated from Equation 
()4.3|) will decrease as /x increases. Another way to reach this conclusion is to think about 
the fi component of the transition charge density. As /x increases, the centrifugal potential 
will keep the transition charge density farther away from the z axis. Thus absorption, which 
occurs close to the z axis, will have less effect on high /x transitions. Since it is absorption 
that gives rise to a minimum effective value of b, less absorption will mean a smaller value 
of b mhl . 

We have done calculations for /x = 2, and found a decrease in fe min , compared to the /x = 1 
values presented in the last section, by about 0.05 fm when the projectile is 208 Pb. For the 
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range of targets and bombarding energies we have studied, this decrease in 6 min is equivalent 
to an increase in CEX cross-section of about 0.7% to 1.5%. This difference is small enough 
to be ignored in practical calculations. If the projectile is 16 0, the corresponding error in 
6 m i n is about 0.1 fm, which is equivalent to errors of about 1.5% to 7% in CEX cross-sections. 
Thus, for lighter projectiles, and in situations in which fi > 1 transitions are expected to be 
important, the values of 6 min calculated with Equations ()4.4|) . ()4.5|) should be interpreted 
only as upper limits. 

To do a calculation that includes the effect of absorption accurately when several values 
of \x are important, it is necessary to give up the picture implied by Equation (jl.2j) . in which 
absorptive processes occur independently of the Coulomb excitation. Rather, the coupled 
equations that are used to calculate the transition amplitude would have to include absorp- 
tive processes along with electromagnetic processes. This would be analagous to the way 
that optical model analyses of inelastic scattering employ an imaginary potential to simulate 
absorption into other channels in the calculation of the inelastic scattering amplitude. 

APPENDIX A: SOME BASIC FORMULAE OF THE SEMI-CLASSICAL 
APPROACH TO RELATIVISTIC COULOMB EXCITATION. 

The Fourier transform of the matrix element for the transfer of angular momenta (A, /i) 
in the target transition (fi a — > <fip, due to the time-dependent electromagnetic field of a spher- 
ically symmetric projectile moving with speed v along a trajectory with impact parameter 
b is 

00 dt' 
2Z P e 



Vp a {u,b) EE / —Vp a (t',b) 



hv 



-e 



Gx,» J dV (p^(r') - ■ j^r')) jx O^A r A 



x (-) 
Here pp a {v') and j/fa(r') are the target charge and current transition charge densities for the 
states (p a — > (f)p, and is a modified Bessel function. The coefficients are defined by 

ft* = ^ (^f) " v / 4 7 r(2A + l)(A- /U )!(A + /U )! 



x 



^ (2 7 ) 2 ™(n + y u)!n!(A-/i-2n)' ^ 
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This expression assumes there is no overlap between projectile and target charge. The time 
reversal phase convention is used, so that the spherical harmonics has an extra factor of 
i compared to a Condon- Shortley spherical harmonic. 

If it is assumed that no contribution to RCE occurs from b < b m \ n and no nuclear inter- 
actions occur for b > b m i n , then the first-order time-dependent perturbation theory approx- 
imation for the cross-section for the RCE population of a state 0^ in an even-even nucleus 
is 

POO 

a = 2n bdb\Vf3 a (uj,b)\ 2 (A3) 

J 6=6min 

The value of u to be used here is | times the excitation energy of This is sometimes 
referred to as the "on-shell" u value. For the calculations in this paper, the excited state 
corresponds to a one-phonon excitation of the GDR, and thus the on-shell value of to is 
given by Equation Because all the b dependence of Vp a (u, b) is contained in the factor 

K ^ ( ^-J, the b integral in Equation (EqnA.3) can be performed exactly, with the help of 



t2 

2 7 S 



(K^(x)) xdx o 



(A4) 
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Figure Captions 

Figure 1. The reactivity of a projectile nucleon is represented by a grey disk (seen 
edgewise) located at rp relative to the projectile center. Points on the disk are labelled by 
s relative to the center of the disk, and by rp relative to the target center. 

Figure 2. Goldhaber- Teller picture of the GDR, with proton and neutron spheres oscillat- 
ing relative to each other. The oscillation variable is rp^, the vector connecting the centers 
of the two spheres. 

Figure 3. The circles correspond to b min values calculated using Equation ([4.3)1 . when the 
bombarding energy is 1 GeV per nucleon. The crosses are the same, but for a bombarding 
energy of 5 GeV per nucleon. The continuous lines were calculated using Equation ()4.4j) . 
with A = 1.3115 fm and n = 1.0509 fm. 
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Bombarding energy 


&ZZ 




°"ZN 


-E>ZN 


(GeV per nucleon) 


(mb) 


(GeV/c)- 2 


(mb) 


(GeV/c)- 2 


1 


47.849 


5.814 


40.221 


4.1 


2 


45.024 


6.349 


42.973 


5.904 


3 


42.493 


6.847 


42.520 


6.645 


4 


41.307 


7.280 


42.124 


7.384 


5 


40.809 


7.737 


42.567 


8.069 



TABLE I: Parameters to determine the spin-averaged nucleon-nucleon scattering amplitudes, in- 
terpolated and extrapolated from figures given by Igo Q|. 





74 Ge target 


202 Hg target 


Projectile 


<tlw (barns) 


O-overlap (barns) 


ctlw (barns) 


Overlap (barns) 


16Q 


1.9 x 10" 2 


-9.3 x 10" 5 


7.3 x 10" 2 


-9.9 x 10~ 5 


208 pk 


1.3 


-3.5 x 10~ 4 


5.8 


-3.5 x 10" 4 



TABLE II: <tlw and cioveriap, defined in Equation (|3.14j) . for 16 and 208 Pb projectiles bombarding 
74 Ge and 202 Hg targets, at 2 GeV per nucleon. 



cm. 




ic-fr 



P N 




Protons 



v 



